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Abstract 

We prove a powerful scaling property for the extremality condition in the 
recently developed variational perturbation theory which converts divergent 
perturbation expansions into exponentially fast convergent ones. The proof 
is given for the energy eigenvalues of an anharmonic oscillator with an arbi- 
trary ^-potential. The scaling property greatly increases the accuracy of the 
results. 
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In a series of recent papers and a textbook it has been demonstrated that divergent pertur- 
bation expansions in quantum mechanics can be turned into exponentially fast convergent 
ones with the help of variational perturbation theory [1-3]. The purpose of this note is to 
show that the calculation can be greatly simplified by observing that the extrema of the 
energy depend universally, i.e., independently of the coupling strength, upon a simple scal- 
ing variable a. Since a is a polynomial of degree p/2 + 1 in the variational parameter, this 
reduces the degree of the polynomials to be extremized by this factor (3 for the quartic 
potential). The the order of the approximation can therefore be raised by a factor p/2 + 1. 
With the exponentially fast convergence, this leads to a great increase in the accuracy of 
the results. 

Consider a general symmetric anharmonic oscillator with a potential 



The standard Rayleigh-Schrodinger perturbation theory yields for each energy level a power- 

p+2 

series expansion in the dimensionless reduced coupling constant g = g/cu 2 : 



The expansion coefficients ef w are rational numbers which can easily be computed by a 
recursion scheme derived a long time ago by Bender and Wu [4] . A direct summation of the 
series (2) is meaningless due to the factorial growth of the coefficients ef w giving rise to a 
vanishing radius of convergence. An approximate evaluation is only possible at very small 
g after truncating the series at an order I m 1/g. 

An exponentially fast convergent evaluation of the series (2) becomes possible for all 
coupling strengths [3] with the help of variational perturbation theory. It is here where the 
restriction to even powers p in the potential (2) is necessary; otherwise a second variational 
parameter is required. 

The variational perturbation expansion is derived as follows. First, the harmonic term 
of the potential is split into an arbitrary harmonic term and a remainder 



V(x) = —x 2 + gx> (co 2 ,g>0) 



p = even. 



(1) 




(2) 
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Rewriting 

V(x) = ^x 2 + V int (x), (4) 

with an interaction 

V iQt (x) = 9(rx 2 + x p ), r = (u 2 - Q 2 )/2g, (5) 

one performs a perturbation expansion in powers of g at a fixed r, 

i 



E N (g,r) = nf:e l (r) (°) . 

1=0 \fi 2 / 



(6) 



The calculation of the reexpansion coefficients e;(r) up to a specific order N does not require 
much additional work since it is easily obtained from the ordinary perturbation series (2) by 



replacing u by ^/Q 2 + 2gr and expanding in powers of g up to the iVth order. This yields 



BW 

3=0 



f (l-^j)/2^ 



{2rn E ^) l - j . (7) 



V l ~ J J 

The truncated power series 

W N (g, n) := E N (g, (u, 2 - n 2 )/2g) (8) 

is certainly independent of Q in the limit k — > oo. At any finite order, however, it does 
depend on Q, the approximation having its fastest speed of convergence where it depends 
least on Q, e.g., at points where dW^/dfl = 0. If we denote the order-dependent optimal 
value of Q by fl^, the quantity W^(g, fl^) is the new approximation to E(g). 

The extremization yields a large number 3N of il-values from which to select the best 
Qn by choosing the flattest extremum, i.e., the one with the smallest second derivative. We 
shall now demonstrate that the scaling variable 

a = -2rQ— = ^ '- 9 
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makes the extremality condition dW^/dQ, = a universal function of a: 



P N ((J) = 0, 



(10) 



where P N is a polynomial of degree N in a. The polynomial is obtained by forming the 
derivative of Wjs;{g,Q) and removing the iVth power of the reduced coupling constant as 
follows: 



dW, 



N 
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N 



(11) 



The condition (10) increases the order to which the expansion can be evaluated by a factor of 
p/2+1. Since the convergence is exponential, this greatly increases the accuracy which can be 
reached by a variational perturbation expansion. This nontrivial property of the variational 
approach was discovered empirically while treating the quartic anharmonic oscillator in 
Ref. [5] (for details see [3]), but not understood. 
We prove Eq. (11) and derive the formula 

de N+1 (a) 



Pn(<t) 



da 



Together with (7), we obtain the polynomial of degree N: 

de N+ i{a) 



Pn(o) 



N 



2£' 

3=0 



da 



,BW 



' (i - N 

[ N + l-j J 



(N + l-j)(-a) 



N-j 



The proof proceeds as follows. Differentiating 



with respect to Q we obtain 



dW, 



N 



dn 



E 

1=0 



p + 2 



de. 



2 te 'W + n 5si 



(12) 



(13) 
(14) 



(15) 



(16) 



With the help of 



n 



de t 



n 2 ? p-2 ' 

2 + a 

9 2 



de t 
da 



(17) 



this becomes 



dW. 



N 



dfl 



N 

£ 

1=0 



2 J lK ' \ g 2 I d<r 



>i2 2 



(18) 



Rearranging the sum gives 



dW N 2 de f g 



dQ ~ do 

N-l 

+ £ 



i=i 



p + 2\ p-2 de, de l+1 



£+2 



(19) 



+ 



, P + 2 \ p-2 de N 
1 ; — N e N H — <r 







2 "V ; 2 " da 

The first term vanishes trivially since e does not depend on a. 

The crucial observation essential to our result is that the sum in the second line vanishes 
term by term: 



p + 2 \ p-2 de t , n de l+1 



(20) 



To derive this we take the derivative of (7), 



de 



i+i 



da 



2£< 

3=0 



,BW 



V 



(1 " V -f3)l2 
l + l-j 



{3-1- 1) (-*) 



(21) 



and use the identity 
( 



(1 " *?3)/2 
1 + 1-3 



\ 



J 



^3 + 21-1 
3-1-1 



( 



(i - *%m 

1-3 



(22) 



to rewrite (21) as 



de, 



i+i 



da 



j=0 



BW 



(1 - P ~f3)/2 
1-3 



\ 



J 



p-2 



j + 21 - 1 (-a) 



(23) 



Since 



■5 



p 



- 2 de t 

a— — 

2 da 



( 



3=0 



BW 



(i - *fm 

1-3 



p 



1-3 



(24) 



we find 



p - 2 dej rfe ;+ i 

cr— h z — - — 

2 acr acr 



p + 2 



I -I )e l . 



(25) 



This proves, indeed, that each term in the sum of (19) vanishes identically. 

Finally, inserting Eq. (20) for / = N into the last term of Eq. (19), we obtain the desired 
result (11): 

dW N 



N 



dn 



9 \ 2 de N+1 (a) 



p+2 



da 



(26) 



The discovery of the scaling variable a was essential for the recent determination of the 
strong-coupling expansions of the energy levels of the quartic anharmonic oscillator to a 
great accuracy [5] (described in detail in [3]). 

It will be interesting to see whether there exists similar universal scaling variables for 
odd interaction powers p where the asymmetry of the potential requires the introduction of 
a second variational parameter, a shift in the average position of the harmonic trial potential 
in [6]. 
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